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Observation of an anomalous temperature dependence of the spin susceptibility, along with a 
spin gap in NMR, in the quantum paramagnetic normal state of Iron Pnictides is a signature of 
an unusual metallic state. We argue that both these anomalous features are associated with a wide 
fluctuational regime dominated by dynamical, short-ranged and frustrated spin correlations in a 
strongly correlated metal. Using LDA+DMFT, we show that both these features can be quantita- 
tively undertstood in the doped Fe-pnictides. We argue that such spin correlations naturally arise 
in a Mottness scenario, where an effective, dualistic description involves coexisting renormalized 
quasiparticles and effectively localized moments, arising from the same set of d-bands. 

PACS numbers: 74.70.-b, 74.25.Ha, 71.27. +a 
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The precise mechanism of unconventional supercon- 
ductivity (U-SC) in the recently discovered Iron Pnictides 
(FePn) is presently a hotly debated issued U-SC arises 
on the border of a geometric frustration (GF)-induced 
"stripe" -spin-density- wave (SSDW). 2 Short-ranged, frus- 
trated spin correlations, expected to survive doping- 
induced destruction of SSDW order j 3 - are then expected 
to play an important role in the normal and SC phases 
of FePn. 

Study of magnetic fluctuations in an U-SC can help un- 
earth the symmetry of the SC order parameter, as shown 
by detailed studies for cuprates^. In FePn, extant exper- 
iments already show anomalous behavior. Specifically, 
NMR studies reveal normal state pseudogap behavior— 
and U-SC, as seen from opening up of a spin gap at 
T* = 150 - 200 K, and Tf 1 ~ T" with n = 2.2 - 2.5 
for T « T c . However, other probes reveal anisotropic, 
albeit fully gapped, structure of the in-plane gap func- 
tion. While this is the expected response of a GF spin 
system in its quantum disordered phased FePn are met- 
als, albeit presumably close to a Mott- Hubbard instabil- 
ity^ Are there other signatures of anomalous magnetic 
correlations in FePn? Extant data reply in the affir- 
mative: the uniform spin susceptibility shows a linear- 
in-T dependence for T > Tsdw(x) up to 800 K for 
both the 1111- and 122-FePn. 6 ' 7 This anomalous behav- 
ior persists upon destruction of the SSDW, even going 
over to a x(T) — T 1+n form at lower T, with n < 1. 
Actually, just above T c for the doped La-based 1111- 
FePn with x > 0.05, the spin susceptibility obeys the 
x{T) = Xo + AT 1+n law, and the linear-in-T dependence 
is recovered smoothly only at higher T. This observa- 
tion is strong circumstantial evidence for relevance of GF, 
since the x(T) — T law is also seen above Tjy in another 
GF system, Nao^Co02, which is also a poor metal like 
doped FePn, 8 However, Cr and its alloys, widely believed 
to be described in a more itinerant framework, also show 
the linear-in-T susceptibility. The distinguishing feature 
of the FePn with x ~ 0.1 is, however, that they exhibit 
the spin gap in NMR in the quantum disordered phase. 
This observation puts them into the "strong coupling" 
category, since the spin gap in the paramagnetic phase 



is beyond reach of an itinerant description. A consis- 
tent theory of magnetic fluctuations in FePn must then 
reconcile the spin-gap in NMR with X.(T) ~ T 1+n , with 
< n < 1 in a single picture. 

Earlier theoretical studies have focussed on both itiner- 
ant weak- and strong coupling routes?^ In the itinerant 
view, a x(T) ~ T requires nesting of the electron- and 
hole-like Fermi sheets. This involves fine tuning, since 
no such nesting can be invoked for doped 1111-FePn, 
even as the x(T) ~ T behavior persists. The second 
view, based on (presumably) Kondo-like coupling of itin- 
erant carriers to effectively localized spins in a J\ — J2- 
Heisenberg model, does not require fine tuning, but is 
strictly valid only in the localized regime. 2 Experimental 
evidence, however, places the FePn in the intermediate- 
to-strong correlation regime^ 2 - where the dualistic aspect 
of self-consistently correlated d electrons should manifest 
itself. In this Mottness scenario, \both the "itinerant" 
carriers and the "local" moments arise from the same set 
of d-electrons, simultaneously giving renormalized LDA- 
like "bands" as seen in dHvA 10 and "localized" magnetic 
responses seen in NMR and INS studies. We thus use 
LDA+DMFT, which is known to capture "Mottness" ex- 
actly at the mean-field leveh 11 ' 12 

In a study of magnetic fluctuations, the central quan- 
tity of interest is the dynamical spin susceptibility, 
x(q, ui) — J2 a b Xa,b(q, w), where a, b are all d-orbital in- 
dices, and q, tu are the momentum and energy transfers 
in INS. Viewing FePn as strongly correlated systems, we 
construct x(q, w) in terms of the full LDA+DMFT prop- 
agators computed in earlier worki 13 ' 14 ' 15 Good quanti- 
tative agreement between LDA+DMFT and key exper- 
iments in both, the normal and U-SC states, has been 
shown there, supporting our choice. The prescription 
is to replace the band Green functions used in weak- 
coupling RPA-like approaches 1 ^ by their LDA+DMFT 
counterparts. This ensures that the dynamical aspect of 
strong, local, multi-orbital (MO) correlations is included 
from the outset. 

For the incoherent "normal" metallic state in FePn, 
after replacing the bare G aa (k,uj) with G aQ (k, u>) = 



QLDA+ 



DMFT 



(k,uj) = [w - e ka - T, a (u>)Y 



and intro- 
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ducing the spin operator <S a , M (q) = 

q)°^, CT ,<7' c a,At,^( k )) with V = X >V> 

cal spin susceptibility reads 



2 ^a,fi,a\ 

the "bare" dynami- 



Xo^.tfa) w ) = K li J dv J de P™°( k + ^ ' 

k,cj ; 

n F {v) - n F (e) 



p bba > (k, e) 



uj + v — e + irj ' 



where K™ 



2 a,aa' 



a are Pauli 
p Q (k, v) = 



<,^>> the 

matrices, and rip the Fermi function. 
(— l/7r)ImG QQ .(k, f) is the one-particle spectral func- 
tion for orbital a. Including the ladder vertex 
in an infinite summation of "ladder" diagrams us- 
ing RPA, the renormalized magnetic susceptibility, 
Xa,b(q,w) = [xo^ b (q,w) - J(q)]~\ where, following^ 
Xo,a,b(q, is evaluated in DMFT. The "bare" bub- 
ble contribution is found as x^(0> w ) — X ( w ) = 
CY.a,b5 dv Pa{v)pb{u + v)[f{v) -f(u + v)], and J(q) = 
Ji(cos(g x a)+cos(gj,a))+ J2COs(g x a)cos(<7j / a), with J\ ~ 

t 2 t' 2 
jji'+j and >^2 — u'+Jh being the frustrated superex- 

change scales in FePn^ Using Xo,o,6(q>w), the dynami- 
cal spin susceptibility x"(Q: can be now expressed in 
terms of the full DMFT propagators computed earlier for 
the incoherent "normal" stated 

Using x(q, u>) at low T, the uniform (q = 0) spin sus- 
ceptibility is now estimated as 



X(T)= X (0,T) 



+ OC 



Notice how frustration appears in the "bare" bubble com- 
prised of fully renormalized DMFT propagators, which 
contain contributions from the (frustrated) hoppings (via 
DMFT propagators), as well as in J(q) via the RPA sum. 
In a companion work^ we have shown how this yields 
nice agreement with the NMR relaxation rate over the 
full T range for FePn,— including spin gap formation and 
unconventional power-law form of T^ 1 below T c . Here, 
we extend that study to investigate the spin susceptibil- 
ity in detail. We explicitly demonstrate that x(T) ~ T 
at high T, with a smooth crossover to a xCO — T 1+n 
form at lower T with n smoothly increasing from zero to 
unity as T is lowered, for the doped FePn which become 
superconducting at lower T. We show how both NMR 
and susceptibility data in the "normal" state can be ra- 
tionalized in terms of strong, frustrated spin correlations 
arising in a multi-band, correlated system on the verge 
of Mottness. 

Before presenting our results, we make a few relevant 
remarks. In contrast to earlier studies;^ our model in- 
corporates all five d bands, and explicitly captures the 
anisotropic three-dimensional band structure of the 1111- 
FePn. This should become relevant at low T, and a 
smooth crossover to D = 2 physics should occur with in- 
creasing T above Tsdw(x)- In such a D = 2 regime, one 



might associate the x(T) — T law with the fluctuational 
regime that occurs below the mean- field (Tmf — J 2 >> 
TsdWjTsc) crossover^ The wide T range over which 
these latter exist points to dominant effects of non-local 
spatial spin correlations up to rather high T. Combin- 
ing this with the GF nature of FePn and observation of 
the spin gap below T* in the 1111-FePn suggests that the 
spatial extent of these spin correlations should be of order 
a Fe — Fe unit cell. These are precisely the correlations 
accessed by our formulation above, though our neglect 
of vertex corrections of the "crossing diagram" type is 
an approximation. We note, however, that a similar ap- 
proximation is shown to reconcile the ARPES lineshapes 
and the neutron scattering results in cuprateSfii This in- 
dicates that such vertex corrections are small, justifying 
our approximation. 
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FIG. 1: (Color online) T-dependence of the uniform spin sus- 
ceptibility, x(T), for the x = 0.1 La-based Fe-pnictide. No- 
tice the smooth evolution of the linear-in-T law at "high"- 
T > 200 K to a T 1 ' 8 law at lower T, in nice agreement with 
experiment.— Also note that the [/-dependence of the sus- 
ceptibility is strong evidence for a sizably correlated multi- 
orbital scenario in FePn. For clarity, both theoretical curves 
are shifted downward to coincide with the experimental curve 
at T c . 

We now present our results. In Fig. [TJ we show 
the T-dependent uniform spin susceptibility, x(T), for 
LaO\- x FeAsF x at a representative electron doping, x — 
0.1, for two values of U corresponding to intermediate- 
to-strong (U = 4.0 eV) and to weak (U = 2.0 eV) cou- 
pling. Since we do not treat the SSDW instability here, 
we focus only on the regime x > 0.04, where we study 
the T-dependent x{T) up to lower T in the "normal" 
state. For x = 0.1, SSDW order is destroyed by doping. 
However, in full agreement with experiment, x{T) con- 
tinues to show the linear-in-T dependence in the range 
200 K< T < 800 K (only shown up to 350 K), even ex- 
hibiting a smooth increase of curvature (~ T 18 ) at even 



3 



lower T in the vicinity of T c . In fact, for x = 0.1, we 
have fitted x(T) up to 500 K with the 2.03 + 0.42T 1 - 8 
law, quite distinct from the linear-in-T behavior seen for 
small x. 7 ) Only at high T is the linear-in-T behavior re- 
covered for x ~ 0.1. The T-dependence of our computed 
susceptibility is very close to this fit to experiment, as 
shown in Fig. [TJ Clearly, x(T) with U = 2.0 eV sub- 
stantially deviates from experiment, while a much bet- 
ter agreement is obtained with U = 4.0 eV, testifying 
to the importance of sizable multi-orbital electronic cor- 
relations in FePn. For sake of discussion, we show our 
earlier theoretical result for the NMR relaxation rate, 
[TiT]- 1 , 14 with U = 4.0 eV in Fig. H The absence of 
a Korringa form ([TiT] -1 =const) at any T is striking. 
Around T* ~ 150 K, clear, gradual opening up of a spin 
gap, again as seen in NMR work, is found. The correla- 
tion between the spin-gap in NMR and x(T) — T is not 
easily obvious: x(T) continues to vary smoothly even as 
[TiT]" 1 shows the spin gap behavior. 




FIG. 2: (Color online) Low T behavior of the NMR [TiT] -1 
for LaOi- x FeAsF x with x = 0.1, 0.2, corresponding to doped 
samples with no SSDW order. The emergence of the spin 
gap around T* ~ 150 — 200 K for the doped FePn is clearly 
resolved, in good agreement with experiment.— 

Such semi-quantitative agreement, even in details, 
with both NMR and spin susceptibility experiments is a 
very gratifying feature of our modelling. Deeper insight 
into these anomalous responses, and of their connection 
to U-SC in doped FePn, is clearly desirable. To this end, 
we start by noticing that x(T) — T is found in the range 
< T < Tmf — O(J) even for a classical, unfrustrated 
lattice in D = 2M- Since the D — 2 Heisenberg model 
cannot have a finite Tv, it follows that this linear be- 
havior must be associated with the fluctuational regime 
below the mean- field crossover, as discussed by Zhang 
et al. in the localized limit of the J\ — Ji models In a 
GF system like FePn, this regime involves short-ranged, 
dynamical spin correlations above Tsdw ■ I n "bad metal- 



lic" FePn, the "Mottness" philosophy allows us to de- 
scribe the low-energy magnetic fluctuations by an effec- 
tive J\ — Ji models In this model, there is a wide T range, 
Tsdw < T < T Sl where translational symmetry is spon- 
taneously broken but spin-rotational SU(2) invariance is 
noti^ Such a transition has interesting consequences: 

(i) it naturally accounts for opening up of a spin gap, 
leading to consistency with NMR data, and, 

(ii) it breaks the symmetry between q = (tt, 0) and 
(0,7r), already above Tsdw- This breaks the lattice 
rotational symmetry (Ref. y) and should lead to the 
tetragonal-orthorhombic (T-O) distortion,- 9 accompa- 
nied by a striped-SDW long range order at Tsdw- In 
fact, in general, T s > T$Dwr^ supporting such a view. 
This could also lead to a nematic instability, a possibility 
that has already been proposed in the FePn context^ 

(iii) Finally, these dynamical short-ranged spin corre- 
lations should survive the doping-induced melting of the 
SSDW. These nearest- and diagonal neighbor correla- 
tions are attractive candidates for the "glue" that leads to 
short-coherence length U-SC in FePn. In fact, the U-SC 
gap function should then have neartest neighbor (n.n) 
and diagonal components: these exactly correspond to 
the ex-s and s^-pair components. Interestingly, rigor- 
ous analysis requires that both co-exist in FePn.?! This is 
precisely the in-plane, nodeless part of the SC gap func- 
tion that we have proposed in Ref. [TH as an instability 
of the incoherent "normal" state in doped FePn. Since 
the pair wavefunction has only n.n and next-n.n compo- 
nents as above, this also affords a way to rationalize the 
short SC coherence length, high upper critical fields, and 
Uemura scaling, all of which are observed in FePni22 

Hence, in our work, the anomalous spin susceptibil- 
ity arises from effectively localized, short-ranged and 
sizably frustrated spin correlations above Tsdw m the 
1111-FePn. That the FePn are bad metals (the drop 
in p{T) at low T in some of the FePn, especially in 
LaO\- x FeAsF x , actually correlates with the opening of 
the spin gap^ and, cannot be taken as evidence of FL 
behavior )£ means that these short-ranged spin correla- 
tions are further damped by the dynamical charge fluc- 
tuations in reality. As a further check of our analysis, 
we notice that, choosing the spin gap, A ~ 150 — 200 K, 
from our NMR result, the dc resistivity is expressed as 
Pdc{T) ~ p + CTe~ A/fc£iT , in good agreement with the 
observations for x = 0.1 in LaOFeAs; 23 interestingly, 
extraction of A from the dc resistivity gives A = 164 K 
for x = 0.1, in good agreement with our estimate from 
NMR. When the spin gap closes, the dc resistivity re- 
covers its characteristic linear-in-T behavior, as indeed 
seen^ This is again reminiscent of the evolution of 
Pdc(T) in underdoped cuprates with increasing hole dop- 
ing^ In a "Mottness" scenario, with relatively stable lo- 
cal moments (this should not be confused with small 
sublattice magnetization in FePn, which can be drasti- 
cally lowered by frustration)^ the effective Heisenberg 
limit should be qualitatively correct below an exchange 
scale renormalized downward by metallicity. We empha- 
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size that this argument does not constitute a strict map- 
ping from our DMFT results to an effective frustrated 
Kondo-Heisenberg model (for a more formal mapping, 
see Ref. [H); however, the latter does allow deeper insight 
into the DMFT results. 

Finally, at lower T for x ~ 0.1, the dimensional 
crossover from D = 2 to D = 3 should occur. This is 
what leads to the deviation of the low T susceptibility 
from the linear-in-T behavior. Our finding of an approx- 
imate x( T ) - Tl ' S law in this regime (in nice agreement 
with the x = 0.1 La-based 1111-FePn) suggests multi- 
orbital effects at work: in particular, we argue that small 
inter-layer couplings (hopping between d z i orbitals on 
neighboring layers, and concomitant superexchange, J±), 
setting in at low temperatures, lead to anisotropic D = 3 
physics, giving an approximate x(T) — T 2 law, smoothly 
evolving from the linear-in-T variation at higher T . 

Our analysis thus leads to a very different picture com- 
pared to the one obtained from an itinerant view. 9 In 
the itinerant view, spin gap generation requires rather in- 
volved extensions of the HF-RPA approach hitherto used 
therein^ and these should bring the theoretical descrip- 
tion closer to the limit considered here. Also, evolution 
of x(T) ~ T to the T 18 dependence with doping (for 
x = 0.1 and beyond see Ref. @) remains to be clarified 
there. In our "strong coupling" picture, this behavior 
arises from the 2D — 3D crossover which should set in 
with increasing x at low T in the quantum paramag- 
netic metal. Using LDA+DMFT, we show how both, the 
X(T) ~ T at high T, as well as its T 2 -like behavior for low 
T in the x > 0.06 regime, can be obtained in a strongly 



correlated picture. Both the linear-in-Tsusceptibility and 
the spin gap are thus associated with the wide ranged 
(in T) fluctuational contributions associated with short- 
ranged, geometrically frustrated spin correlations in the 
quantum disordered regime of an effective, doped J± — J2 
antiferromagnet. Our view is also fully consistent with 
ARPES, 25 optical^ and /LtSRSi data, all of which are 
consistent with a strongly correlated metal giving way to 
U-SC. 

In conclusion, we have studied the temperature de- 
pendence of the uniform spin susceptibility, x{T), for 
the 1111-FePn using the first-principles LDA+DMFT 
method. Good semiquantitative agreement with pub- 
lished experimental data testifies to the relevance of 
short-ranged and dynamical multi-band electronic cor- 
relations in FePn. The same approach also yields a satis- 
fying description of the NMR [TiT]^ 1 as a function of T 
in an earlier work.— We thus identify both, the spin gap 
in NMR and T dependence of x{T), as manifestations 
of the wide (in T) fluctuational regime of a sizably frus- 
trated magnet between its actual ordering scale (Tsdw) 
and the mean-field crossover (Tmf)- Our study pinpoints 
the crucial roles of strong multi-orbital electronic corre- 
lations and sizable geometric frustration to understand 
magnetic fluctuations in the incoherent metallic state of 
the 1111-Iron Pnictides. 
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